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Abstract 


Sink strength is a fundamental quantity in modeling the microstructure 
evolution of irradiated materials by the mean-field approaches. The analyti- 
cal expressions for different sinks have been extensively studied. The Object 
Kinetic Monte Carlo (OKMC) simulations were subsequently performed to 
corroborate the expressions and to guide the development of analytical the- 
ories. Although a general agreement was found between the theory and 
simulation, there are still some discrepancies in the case of the small spher- 
ical sinks and the dislocation lines. In this work, OKMC simulations were 
performed to study the sink strength of spherical sinks and dislocation lines. 
Our results revealed the origins of discrepancies between the theory and simu- 
lation for small sinks, high volume fraction sinks and periodically distributed 
dislocation lines. The theoretical corrections were also proposed according- 
ly. These corrections can extend the capabilities of mean-field approaches to 
properly reproduce the defect evolution of irradiated materials, and provide 
a better insight into the experimental phenomena, such as the defect cluster 


nucleation, the irradiation-induced swelling and the blistering. 
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1. Introduction 


The degradation of material properties under irradiation is a main con- 
cern in reactor design. Under the high-energy particle irradiation, numerous 
mobile defects such as interstitials and vacancies are created. The clustering 
and aggregation of these defect are generally relevant to the the formation 
of blisters and voids, which leads to the degradation of materials. These mo- 
bile defects can interact with each other and with other sinks such as grain 
boundaries and dislocation lines. The ability of these sinks to capture the 
migrating defects can be quantitatively described by the sink strength that 
was applied in mean-field approaches for simulations of defect evolution. 

The sink strength is proportional to the inverse square of average free 
migration distance of defects. It is connected to the size, shape and density of 
sinks, as well as the dimensionality of defect migration, i.e., three-dimensional 
(3D), one-dimensional (1D), or mixed 3D/1D. Based on the diffusion theory, 
many theoretical work has been carried out to obtain the analytical sink 
strength for different sinks. Brailsford and Bullough et al. have investigated 
the expressions for the spherical sinks [1] and grain boundaries [2] under the 
3D limit. Similar expression for dislocations was derived by Wiedersich [3], 
corrected by Nichols [4], and shows a good agreement with the phase-field 
results [5]. Theoretical sink strengths for these sinks under the 1D limit 
were systematically investigated by Barashev et al. [6]. Recently, the Object 
Kinetic Monte Carlo (OKMC) method [7-9] which is expected to provide the 
sink strengths spontaneously, has been used to corroborate these expressions 
and to guide the development of the theories. The OKMC simulations by 
Malerba and Jansson et al. [10, 11] showed good agreements with analytical 
theory in both the 3D and 1D limits, as well as in the transition regime of 


3D to 1D diffusion. The transition region is described by a single-variable 
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function (master-curve), which was proposed by Trinkaus and Heinisch et al. 
and also confirmed by the OKMC simulations [12-15]. 

Despite of the general agreements found between theoretical works and 
simulations, unexplained discrepancies under certain situations have also 
been found in previous studies. More specifically, the OKMC simulation 
results always show an underestimate of the analytical sink strength for s- 
mall sinks [10, 13, 16], and a huge overestimate for high volume fraction sinks 
(10, 11, 13]. An irregular discrepancy for periodically distributed dislocation 
lines was also pointed out[11]. Heinisch et al. [13] attributed the underesti- 
mation to different nature between the diffusion mechanisms in OKMC and 
theoretical model. They thereby provided a fitted correction term, which 
works well for the sink radii above 2.5a9 (ao is the lattice constant). How- 
ever, they did not consider the effectiveness of the correction for the smaller 
sink (< 2.5a9, such as vacancy clusters containing about 100 vacancies). 
More importantly, the corrected expression cannot describe point defects 
like mono-vacancy because it predicts a zero sink strength. Since OKMC 
is a priori expected to provide the same sink strength to that provided by 
the diffusion theory, additional corrections should be required for OKMC or 
mean-field simulations of defect cluster nucleation. 

In this work, we calculated sink strengths of different sinks through both 
the theoretical and OKMC approaches, for the defects characterized by a 
varying motion dimensionality from the fully 3D limit to the pure 1D lim- 
it. Particular cares were taken for the small sinks, the high volume fraction 
sinks and the dislocation lines in the 1D migration limit. The theoretical 
corrections thereof were also provided to eliminate the discrepancies men- 
tioned above. It was found that much better agreements between theory and 


simulation can be achieved after the corrections. 
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2. Computation method 


An OKMC program is designed to calculate sink strengths of spherical 
sinks and dislocation lines. The entire procedure is similar to that used in 
refs. [10, 11]. Briefly, defects and sinks are treated as objects with specific 
positions in a simulation box and with associated reaction volumes. Proba- 
bility of a thermally activated event, such as defect migration, is given in the 
form of Arrhenius frequency, [; = viewp(—Eai/keT). Here v; is the attempt 
frequency for event i, Ea; is the corresponding activation energy, kg is the 
Boltzmann constant, and T is the absolute temperature. An event is random- 
ly chosen according to its probability. The binary search algorithm is applied 
here to accelerate the choosing procedure. After the event is executed, the 
object configuration will be updated, and the simulated time will be increased 
according to the resident time algorithm [17], Ar = 1/ Oe ae), 
where Nin: and Next are the number of internal and external events, respec- 
tively. Here, two internal events are considered, including defect migration 
and rotation of migration direction. No external events are considered in this 
work, i.e., Next = 0. Note that, when two objects overlap geometrically, reac- 
tions between these two defects will take place. For instance, absorptions of 
a point defect take place when its distance to a sink is smaller than the sum 
of their capture radii. This kind of event will bring no time increment to the 
system. During the simulation, events are chosen and executed repeatedly 
until certain time or steps are achieved. In this work, in order to achieve 
a statistical convergence, simulations will not be stopped until 10000 point 
defects are absorbed. 

Similar to previous works/10, 13, 14], a model handles only idealised sit- 
uations is implemented. In this model, a point defect with zero radius is 


introduced at a random position and allowed to migrate. After its absorp- 
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tion, a new point defect will be introduced instantaneously so that only one 
point defect at a time is presented in the system. It should be pointed out 
that defects created inside spherical sinks or dislocation lines are not taken 
into account. In this way, the sink strength, k?, can be given by: 


2n 
k? = —— 1 
d <m >’ (1) 


where n is the dimensionality of point defect migration, which is taken as 3 
even for 1D migrating in order to consistently trace the transition between 
1D and 3D migration. d is the migration distance of a point defect, which 
equals V/3a9/2 in bec tungsten with lattice constant aj = 0.3165 nm. < m > 
is the average number of jumps before it is absorbed by sinks. Similar to 
ref. [10], the dimensionality of defect motion is changed by assigning a ro- 
tation energy E, to defects, whereby the probability of changing migration 
direction is expressed as exp(—E,./kpT). One can set E, = 0 to acquire pure 
3D migration or set E, to a relatively high value to get pure 1D migration 
(<111> directions in tungsten). In this work, at the chosen simulation tem- 
perature of 723 K, E, = 1.5 eV is enough to provide a fully 1D path, while 
the intermediate values are considered to give a mixed 1D/3D migration. 

A non-cubic periodic box with a fixed size ratio lẹ : ly : lz = a: V2 : 
1 is employed in order to correctively describe the 1D migrating of point 
defects [10]. The irrational ratio of box sides was used to ensure that any 
super-box of this box would be non-cubic. Spherical sinks are randomly 
introduced in this non-cubic box with l, = 350 nm. The radius of spherical 
sinks, R., varied from 0.27 to 10.25 nm. And their number density, N,, 
varied from 1075 to 1.5 x 1074 nm~? for the 3D and 1D migration. For the 
mixed 1D/3D migration, a fixed R, = 4.25 nm was used and N, varied from 
3 x 1075 to 1.5 x 10-4 nm™’. As shown in Fig. 1, the dislocation line is 


introduced as a cylindrical sink whose two opposite faces touch the faces of 
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the simulation box. Only one dislocation line is introduced along l, direction 
of the simulation box. However, under the periodic boundary condition, this 
model is practically equivalent to a regular array of infinitely long dislocation 
lines. Different surface densities can be obtained by changing the box size. 
Here, l, varied from 15 to 45 nm while keep the size ratio fixed, corresponding 
to surface densities from 5 x 1074 to 4.5 x 1073 nm~?. The dislocation radius 
used here, Rg, varied from 0.25 to 10 nm. Some large Ra are not used for 
small boxes due to the box size limit, and some small ones are not used for 


large boxes owing to the computational time limit. 


3. Results and discussion 


3.1. Spherical sinks 


3.1.1. Results without any correction 
The analytical sink strength of spherical sinks under 3D migration limit 
is given by [1]: 
k3 s = 40R,N,(1 + kg,sRs). (2) 


This equation is usually solved by the iterative method. The n‘? order of 
iterative result is calculated by KB sn = ATR: N:(1+ k3 sn-1Rs), with the first 
order given as ksi = 47R,N,. As proposed by previous work [10], it is 
necessary to use high orders of iterative approximation for a good agreement 
between the theoretical expression and the OKMC simulation. For Eq. (2), 
the larger the R, is, the poorer the convergence will be. Even so, as shown in 
Fig. 2, for the largest R, = 10.25 nm used in this work, a good convergence 
is already achieved at the third order approximation. In this work, the fifth 


order approximation is adopted for all calculations using Eq. (2), and for 


similar equation Eq. (4). In the case of 1D migrating defects, expression of 
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the sink strength is given by [6]: 
kis = 6(T RNa}. (3) 


In our OKMC simulations, spherical sinks were introduced randomly but 
without overlap with each other. Under the 3D migration limit, the radius 
of spherical sink, R,, varied from 0.27 to 10.25 nm (0.27, 0.75, 1.25, 2.25, 
4.25, 6.25, 8.25, and 10.25 nm), and their number density, N;, varied from 
3x 10-5 to 15 x 10-5nm~? (1.0, 3.0, 6.0, 9.0, 1.2, and 15.0 x 10-°nm73). The 
two highest densities were not considered when R, = 10.25 nm, because of 
the limitation of non-overlap. For the 1D migration limit, it is impossible to 
produce all the data points explored in the 3D case owing to computational 
time reasons. The same densities to the 3D case are all taken into account 
for R, > 4.25 nm. While only one number density (N, = 107? nm~) is 
considered when R, < 4.25 nm (0.75, 1.25, 1.75, 2.25, 2.75, 3.25 and 3.75 
nm). 

Fig. 3 presents results of sink strength under the 3D and 1D migration 
limits. Similar to previous studies [10, 13], a general agreement is found 
between analytical theory and simulation. Whereas two significant discrep- 
ancies can also be identified for small sinks under the 3D limit and for the 
high volume fraction sinks. These discrepancies can be clearly displayed by 
relative errors, defined as the percentage ratio of simulated-analytical sink 
strength difference to analytical sink strength, i.e., (k2m —Keng)/K2ng X 100%. 
As shown in Fig. 4(a), the simulated sink strength is about 33% lower than 
the analytical one for the smallest sink under the 3D migration limit. Mean- 
while, errors for the same sink radius are close. This indicates that the sink 
radius is the key variable of such errors. However, for the 1D limit that 
shown in Fig. 4(b), similar errors for small sinks are not observed. With 


the growing of sink radius, in the volume fraction range of 1074 ~ 1071, the 
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relative error reduces to nearly zero in both cases, suggesting a good agree- 
ment between theoretical works and simulations. When the volume fraction 
is above 1071, the relative error grows very rapidly with the increase of sink 
volume fraction. This implies a large discrepancy between theory and sim- 
ulation data. In order to eliminate the discrepancies, further investigations 


and additional corrections are required. 


3.1.2. Correction for small sinks 

The discrepancy for small sinks under the 3D migration limit is often 
attributed to different nature of diffusion mechanisms between the theory and 
the OKMC simulation [13, 16]. As shown in Fig. 5(a). For the theoretical 
model, the diffusion is a continuous process, where point defects are absorbed 
exactly at the sink surface. For the OKMC simulation, the diffusion is a 
discrete process composed of a series of individual hops, where the absorption 
takes place when a defect jumps into the sink. Since the jump distance d 
is not zero, point defects in OKMC simulation can penetrate into the sink. 
The penetration area is not an ideal spherical surface but a spherical shell 
(gray area in Fig. 5(a)). Therefore, R, in the Eq. (2) should be revised by 


adding a correction term AR: 
kj, = 4n (R, + AR)N,[1 + kz, (Rs + AR). (4) 


Here, Rs + AR, is denoted by Rey, namely the effective radius of the sink. 

For the 3D limit, the thickness of the penetration area in the OKMC 
model is the jump distance d. Since the sink strength of spherical sinks varies 
linearly with the radius for small sinks, the correction term AR should be 


weighted average of the penetration depth p: 


AR = -p =r — R,, (5) 
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where r = R, — p is the distance to sink center. For a certain r, its possible 
start point (red dot) of the last jump falls on a spherical cap surface with 
radius d centered at destination point (black dot). Form r to r+dr, the pos- 
sible start point is given by shifting the spherical cap surface by dr (shadow 
area in Fig. 5(a)). The volume of this area is given by the differential of 
spherical cap volume|18}: 


di! + (r? — R2)? — 2d2(r? + R2) 
T 


dV = T 


dr. (6) 


The correction term is given by: 


Rs — 2dV 
ApS ee ralf — Re)Anr gdr 


Rs dV 
24V 
Toog Arr? dr 


(7) 


where 4nr? “dr represents the set of start points that can jump to a spherical 

shell between r and r + dr, namely, the weight of r. Substituting Eq. (6) 

into Eq. (7), we have: 

B 5d? + 16dR, — 30 R? 
16d? + 30d R, — 80R? 


AR = d. (8) 


By substituting this correction term into Eq. (4), we calculated percentage 
errors after the effective radius correction. The modified sink strength and 
its corresponding relative error are calculated and summarized in Fig. 6. It 
clearly shows that that discrepancies for small sinks are perfectly removed 
after the effective radius correction. 

We also noticed that there is no discrepancy for small sinks under 1D 
limit. This can be understood as follows. As shown in Fig. 5(b), under 1D 
migration limit, the sink strength is connected to the biggest cross-section 
of the sink perpendicular to the migration direction. Under the continuous 
diffusion mechanism, point defects are absorbed exactly at the sink surface. 


Therefore, the radius of the biggest cross-section is equal to the sink radius 
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R,. While in the OKMC simulation, as stated before, point defects can 
penetrate into the sink. The penetration area here is the gray area between 
the two spherical surfaces separated by d in the migration direction. Res 
is thereby given by the radius of the biggest cross-section of the penetration 
area. Since defects penetrate along the migration direction, Reff is also 
equal to R,. In this way, the discrete diffusion mechanism will not change 


the effective radius Ref, and the correction term AR = 0. 


3.1.8. Correction for high volume fraction sinks 

The discrepancy at the high sink volume fraction region is likely a re- 
sult of different sink configurations. In the OKMC simulation, sinks are 
randomly distributed in the simulation box and do not overlap with each 
other(non-overlap configuration). While in the theoretical model with the 
mean-field approximation, sinks are considered to be randomly and uniform- 
ly distributed in the system without taking care to avoid the overlapping 
between sinks (overlap configuration). In order to verify this speculation, 
we also constructed an overlap configuration in the OKMC simulation. The 
sink strengths under this configuration are calculated and their correspond- 
ing relative error is presented in Fig. 7. Compared with Fig. 4, the relative 
error becomes nearly zero throughout the considered range, which supports 
our hypothesis about the sink configurations. 

The distinct behaviors between two configurations are likely to originate 
from the sink volume fractions. For the non-overlapping configuration, its 


volume fraction f, is given by: 


Due to the overlap between sinks, this formula does not apply to the overlap 


configuration. In the overlap configuration, a larger sink density Ns jap > Ns 
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is required to occupy the same f, to the non-overlap configuration. In this 
case, we define a dummy volume fraction that given by simply summing all 


sinks without caring about the overlap: 
Fy = 40 R$ No,tap/3. (10) 
The relationship between fy and F, is given by differential equations: 
fof + dh) = folFo) + [1 - fil) dF, (11) 


F, (fo =0)=0, (12) 
where 1 — f,(F,) represents the volume fraction of unoccupied regions. So- 
lution of these two equations is: 


1 
1— fy 


With this relationship, we find that the error function fits best with the 


(13) 


F,=I1n 


relative error at the high volume fraction of Fig. 4(b) and Fig. 6(b) is given 


- ‘(inp ff 
e= Bt = (14) 
is i; 


In the low volume fraction region, the probability of overlap is very low, 


and F, is practically equal to f,. Therefore, the relative error arising from the 
sink distribution is invisible. However, as the volume fraction increases, the 
overlap probability grows very rapidly. Therefore, the relative error becomes 
significantly large at the high volume fraction region. It should be pointed 
out that the overlap configuration is unrealistic in real materials, because F, 
can exceed 100%. In order to handle the defect interaction under the high 
sink volume fractions, according to Eq. (14), the theoretical sink strengths 
need to be corrected by multiplying a coefficient 


In = 
( TA J4, 
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In addition, in the previous work of sink strength theory|1], Brailsford et 
al. pointed out an intrinsic part of their method is that the sample volume 
must be large compared with the inter-sink spacing in random arrays. This 
implies that the theoretical work is premised on the assumption of low sink 
volume fraction. Therefore, the additional correction for the high volume 


fraction sink should be reasonable. 


3.1.4. Transition from 3D to 1D regime 
The sink strength of spherical sinks for the mixed 1D/3D migration, 
ki 3.5) can be described by a single-variable function [12] interpolating be- 


tween the 1D and 3D limiting cases, i.e., 


y= s(t Vitae), (15) 


where y and x are two dimensionless variables. In the OKMC simulation, y 


and x are defined as: 


y= ER N (16) 
2 4:2 4 
= lon kis ki (17) 
12° ey 


where len = d\/exp(E;./kgT) is the average migration length between two 
direction changes. 

In Fig. 8 (a), the simulated sink strengths of different densities of sinks un- 
der the mixed 1D/3D limit are plotted as a function of rotation the energy F,. 
All simulations are conducted under the non-overlap configuration. Because 
we do not know the effective radius correction AR under this situation, the 
radius of sinks is chosen to be 4.25 nm, which is large enough to avoid the 
error introduced by effective radius change and also small enough to avoid 
the error from the overlap configuration. Sink density N, varies from 3 x 1075 


to 1.5 x 10-4 nm~?. As we can see, the sink strength smoothly transits from 
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3D to 1D migrating with the increasing E,. And a rotation energy of 1.5 
eV is high enough to ensure a pure 1D migration. The master-curve defined 
by Eq. (15) is illustrated in Fig. 8(b). Here we find that the analytical 


master-curve is in a perfect agreement with our simulated data. 


3.2. Dislocation lines 
According to the previous study [11], the theoretical sink strength of the 
3D and 1D migration defects being absorbed by dislocation lines are [3, 4, 6]: 


2 2mpa(l — meaR3) 
k3 4 = 2 2 , (18) 
0.257 pa R3(4 — Tpa R4) — 0.75 — In (Ray/TPpa) 
and 
k? q = 6 x (Rapa), (19) 


respectively. Here, Ra represents capture radius of the dislocation line, pa = 


I>? is its surface density. Its volume fraction is given by f, = 7R2pgq. 


3.2.1. 8D limit 

Fig. 9(a) shows the sink strengths of dislocation lines under 3D migration 
limit. Similar to the situation of spherical sinks, a good agreement between 
simulation and theory is found but two discrepancies appear for dislocations 
with small Rg and with high volume fraction. The former discrepancy arises 
from the difference between diffusion mechanisms, which are continuous in 
the theory while discrete in the OKMC. This discrepancy can be well recon- 
ciled by introducing the effective radius correction to Eq. (18). The corrected 


equation is given by: 


= 2 pall — mpaq(Ra + AR)?] 
34 Inpa(Ra + AR)?[4 — wpa(Ra + AR)?] — $ — In [(Ra + AR) yraa] 


(20) 


Similar to Eq. (2), Eq. (18) is approximately a linear function of Ry when 


Ra is small. Therefore, the correction term AR should be similar to that 
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for spherical sinks. A detailed derivation of AR is a little bit tedious. For 
simplicity, we shall use the AR from Eq. (8)(simply replace R, with Ra). 
As shown in Fig. 9(b), after the effective radius modification, the relative 
error in the low volume fraction region becomes nearly zero. On the other 
hand, the discrepancy at the high volume fraction cannot be explained as we 
discussed for the spherical sinks. This might be because that the analytical 
expression is based on an assumption that pg is a small value, i.e., the low 
sink volume fraction. And the periodic boundary condition, which might 
cause an interplay between vicinal dislocations, is obviously not included in 
the theoretical model. Since such discrepancy is not obvious until the volume 
fraction is higher than 0.1, this is acceptable because the typical value in real 


materials is usually much smaller than 0.1. 


3.2.2. 1D limit and the influence of the periodic configuration 

Fig. 10 (a) and (b) show the sink strength and the relative error for 
dislocation lines under the 1D limit, grouped by dislocation surface densities. 
No effective radius correction is required here owing to the similar reason 
shown in the Fig. 5(b). The agreement between simulation and theory is less 
satisfactory than the 3D limit. Similar to the case of spherical sinks, a rapid 
growth of the relative error is also observed in the high volume fraction region, 
which can be described by Eq. (14). This indicates that the theoretical 
expression is also based on the overlap configuration. This significant error at 
the high volume fraction region can be eliminated by multiplying a coefficient 
(CFE) to the right side of Eq. (19). 

Figure. 10(c) shows the relative error after the correction. All data 
are positive and show a zigzag behavior as a function of the sink volume 
fraction. These positive errors may origin from a truncation of the tail of 


very long distance migration before absorption. Namely, in a completely 
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random distribution of dislocations, free migration distance of point defects 


before its absorption should follow an exponential distribution [6]: 
P(L) = tRapaexp(—7RapaLl), (21) 


where L is the free migration distance and P(L) is its probability density. 
Any positive L can be achieved with a positive probability density in such 
configuration. However, in a periodic dislocation configuration, such very 
long distance migrations are cut out. Fig. 11 shows a sketch of this cut- 
tail effect. In a periodically extended simulation box, 1D migration defects 
cannot travel to a distance longer than Lew1. Meanwhile, we also find that the 
cutoff distance Leut does not vary continuously with Rg. A slight reduction of 
the dislocation radius from Rg, to Raz might cause a mutation of the cutoff 
distance (Leuti to Leut2). Such cut-tail effect will reduce the average migration 
distance, which consequently causes an overestimate of the theoretical sink 
strength. Therefore, the mutation of the cutoff distance is responsible for the 
zigzag behavior of relative errors. 

Figure 12(a) shows the cutoff distance in OKMC, Leut, as a function of 
the sink volume fraction. The average migration distance, L = (tRapa)~' [6], 
calculated from Eq. (21) is also plotted for comparisons. With an increase of 
the volume fraction, the cutoff distance decreases step by step like a staircase. 
The ratio of Leu to L is shown in Fig. 12(b). Four high peaks and four low 
peaks can be clearly distinguished here. Each of them represents a mutation 
of the cutoff distance. The high peaks in Fig. 12(b) correspond to four 
main peaks of relative error in Fig. 10(c) (labeled with green vertical lines), 
while the low peaks correspond to four shoulder peaks of relative error in 
Fig. 10(c) (blue vertical lines). Although we are aware of the origin of these 
error peaks, it is still very difficult to build a complete random configuration 


to eliminate them. More studies are required in the future to achieve a 
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consistency between theory and OKMC simulation. 

Note that, our results are different from the previous work by Jansson 
et al. [11], where the relative error tends to be negative in the low volume 
fractions, but positive in the high volume fraction region. This difference 
is likely contributed to the different simulation box we used. In our work, 
the ratio of box sides perpendicular to the dislocation is an irrational num- 
ber, ls : ly = V2. While in the previous simulation, the ratios are rational 
numbers. It has been pointed out [10] that in order to prevent endless mi- 
grations of 1D point defects, a non-cubic (non-square in this situation) box 
is instrumental to simulate correctly the 1D migrating defects. However, a 
box with a rational ratio, for instance lẹ : ly = 400a9 : 350aọ, is equally a 
square super-box with Lẹ : Ly = Tle : 8ly = 2800a9 : 2800a 9. The square 
super-box cannot completely prevent the endless migrations from happening 
in low sink volume fraction region. Once an endless migration is started, 
the point defect cannot be absorbed until the simulation is halted manually. 
This should be restrictively forbidden because it would increase the free mi- 
gration distance unpredictively, and lead to a negative relative error in the 


low volume fraction region. 


3.2.3. Transition from 3D to 1D regime 

The transition between 3D migrating and 1D migrating is illustrated in 
Fig. 13(a). The capture radii were carefully chosen to ensure the volume 
fraction f, = 0.005 and f, = 0.03. Because the simulation results fit very 
well with the theory for both 3D and 1D migrating at these volume fractions. 
The transition is more abrupt compared to the simulation of spherical sinks. 
The 1D regime is already reached at E, = 1.0 eV. Meanwhile, in the master- 
curve representation (Fig. 13(b)), the simulation always gives a lower value 


than the analytical one in the 3D migrating region. This was also found 
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in the previous study [11]. However, unlike their results, we find that the 
simulation data of the same volume fraction collapse onto the same curve, 
which indicates that the volume fraction should be the key variable of this 


discrepancy. 


4. Conclusion 


We performed a series of OKMC simulations to study the sink strengths 
of spherical sinks and dislocation lines with different radii and volume frac- 
tions. The simulation results were compared with the theoretical expressions 
derived from diffusion theory to evaluate the validity of these expressions, 
and further, to modify these expressions. 

For the small sinks, the discrepancy between the diffusion theory and 
simulation is caused by the different migration distances used in these two 
methods. This can be resolved by adding a correction term to the effective 
radius of sinks. For the high volume fraction sinks, the discrepancy is at- 
tributed to different sink configurations in two methods. It can be eliminated 
by introducing a compensation term of sink volume fraction to the theoretical 
expressions. For dislocation lines under the 1D migration limit, a zigzag like 
discrepancy is found. It is caused by the periodic distribution of dislocation 
lines adopted in OKMC. 

Based on these results, the origins of the above discrepancies have been re- 
vealed and successfully resolved. This can extend the capability of mean-field 
approaches to properly reproduce the defect evolution of irradiated materi- 
als, and provide a better insight into experimental phenomena, such as the 


defect cluster, the irradiation-induced swelling and the blistering. 
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Figure caption: 

Figure 1 Schematic of the non-cubic box and the dislocation line model in 
OKMC. The periodic boundary condition is applied in all three dimensions. 

Figure 2(color online) The first to the fifth iterative order approximation 
of Eq. (2). Here, the largest sink radius of our simulations, R, = 10.25 nm, 
is used for convergence test. The sink volume fraction is given by f, = 


4r R3N,/3, where N, is the number density of spherical sinks. 
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Figure 3(color online) Sink strength of spherical sinks with different radii 
versus sink volume fraction, lines are results of Eq. (2) and Eq. (3), dots 
are the corresponding simulated results. (a) and (b) represents 3D migration 
and 1D migration limit, respectively. 

Figure 4(color online) The relative error between theoretical and sim- 
ulated sink strength that presented in Fig. 3. (a) and (b) represents 3D 
migration and 1D migration limit, respectively. The solid line represents the 
error estimated by Eq. (14). 

Figure 5(color online) Schematics of a point defect being absorbed by 
a spherical sink under two different diffusion mechanisms. The gray area 
represents the penetration area in OKMC simulation. (a) 3D migration limit. 
The effective radius, Rep, is the weighted average of the penetration depth p. 
Shadow area represents possible start points those can jump into the segment 
from r to r+dr ; (b) 1D migration limit. Ress is the radius of the biggest 
cross-section perpendicular to the migration direction of the penetration area, 
which equals sink radius R, in this case. 

Figure 6(color online) Results of 3D migration defects absorbed by spher- 
ical sinks with different radii, theoretical values are calculated by Eq. (4) with 
AR from Eq. (8). (a) sink strength from simulation (dots) and the theo- 
retical expression (lines). (b) the relative error between simulation and the 
theoretical expression, solid line represents the error estimated by Eq. (14). 

Figure 7(color online) The relative error between theoretical (by Eq. 
(4) with AR from Eq. (8), and by Eq. (3)) and simulated sink strength 
of spherical sinks, under a configuration that sinks are allowed to overlap 
without any reaction takes place. (a) and (b) represents 3D migration and 


1D migration limits, respectively. 
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Figure 8(color online) The transition from 3D to 1D migration limit 
for 4.25 nm spherical sinks, grouped by sink densities. (a) simulated sink 
strength as a function of rotation energy, with temperature T = 723 K. 
(b) the master-curve representation, the solid line represents the theoretical 
relationship of x and y in Eq. (15), dots are the simulation results elaborated 
according to Eqs. (16) and (17). 

Figure 9(color online) Results of 3D migration defects absorbed by dis- 
location lines, grouped by box sizes with a fixed ratio. (a) sink strength from 
simulation (empty dots) and theoretical expression (dotted and solid lines 
are the results from Eqs. (18) and (20), respectively). (b) the relative error 
between simulation and theoretical expression (empty and solid dots are the 
results using Eqs. (18) and (20), respectively). 

Figure 10(color online) Results of 1D migration defects absorbed by 
dislocation lines, grouped by box sizes with a fixed ratio. (a) sink strength 
from simulation (empty dots) and theoretical expression (solid lines). (b) the 
relative error between simulation and theoretical expression, the solid line 
represents the error estimated by Eq. (14). (c) data obtained by removing 
the error described by Eq. (14) from (b). Vertical lines are the corresponding 
position of mutation peaks in Fig. 12(b). 

Figure 11(color online) Schematics of 1D migration in a periodically 
extended simulation box in OKMC, picture is taken perpendicular to the 
direction of dislocation lines. Two kinds of round circles represent dislocation 
lines with different radii. The furthest migration distance, namely the cutoff 
distance, is labeled by Lewy and Lew. Dislocations with radius Ra, (left) 
corresponding to a cutoff distance of Leut1, while a slightly reduction of the 
dislocation radius from Rg to Raz (right) will cause a mutation of the cutoff 


distance from Leys to Leu. 
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Figure 12(color online) (a) the mutation of the cutoff distance (in the 
unit of box size l,), Leut, of 1D migration as the function of volume fraction, 
compared with the average migration distance, L, obtained from Eq. (21). 
(b) the ratio of Leut to L, each peak represents a mutation of the cutoff 
distance. Four high peaks corresponding to the four main peaks of relative 
error in Fig. 10(c), while the four lower peaks corresponding to four shoulder 
peaks of relative error in Fig. 10(c). 

Figure 13(color online) The transition from 3D to 1D migrating sink 
strength of point defects absorbed by dislocation lines, grouped by box sizes, 
where the volume fraction of sinks is fixed at f, = 0.005 and f, = 0.03 
(denoted by empty and solid dots, respectively). (a) simulated sink strength 
as a function of rotation energy, with temperature T = 723 K. (b) the master- 
curve representation, solid line represents the theoretical relationship of x and 
y in Eq. (15), dots are the simulation results elaborated according to Eqs. 


(16) and (17). 
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Figure 1: Schematic of the non-cubic box and the dislocation line model in OKMC. The 


periodic boundary condition is applied in all three dimensions. 
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Figure 2: (color online) The first to the fifth iterative order approximation for sink strength 
from Eq. (2). Here, the largest sink radius of our simulations, Rs = 10.25 nm, is used for 
convergence test. The sink volume fraction is given by fy = 4T RÈN, /3, where N, is the 


number density of spherical sinks. 
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Figure 3: (color online) Sink strength of spherical sinks with different radii versus sink vol- 
ume fraction, lines are results of Eq. (2) and Eq. (3), dots are the corresponding simulated 


results. (a) and (b) represents 3D migration and 1D migration limit, respectively. 
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Figure 4: (color online) The relative error between theoretical and simulated sink strength 
that presented in Fig. 3. (a) and (b) represents 3D migration and 1D migration limit, 
respectively. The solid line represents the error estimated by Eq. (14). 
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Figure 5: (color online) Schematics of a point defect being absorbed by a spherical sink 
under two different diffusion mechanisms. The gray area represents the penetration area 
in OKMC simulation. (a) 3D migration limit. The effective radius, Reff, is the weighted 
average of the penetration depth p. Shadow area represents possible start points those 
can jump into the segment from r to r+dr ; (b) 1D migration limit. Resp is the radius of 
the biggest cross-section perpendicular to the migration direction of the penetration area, 


which equals sink radius R, in this case. 
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Figure 6: (color online) Results of 3D migration defects absorbed by spherical sinks with 
different radii, theoretical values are calculated by Eq. (4) with AR from Eq. (8). (a) sink 
strength from simulation (dots) and the theoretical expression (lines). (b) the relative 
error between simulation and the theoretical expression, solid line represents the error 


estimated by Eq. (14). 
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Figure 7: (color online) The relative error between theoretical (by Eq. (4) with AR 
from Eq. (8), and by Eq. (3)) and simulated sink strength of spherical sinks, under a 
configuration that sinks are allowed to overlap without any reaction takes place. (a) and 


b) represents 3D migration and 1D migration limits, respectively. 
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Figure 8: (color online) The transition from 3D to 1D migration limit for 4.25 nm spherical 
sinks, grouped by sink densities. (a) simulated sink strength as a function of rotation 
energy, with temperature T = 723 K. (b) the master-curve representation, the solid line 
represents the theoretical relationship of x and y in Eq. (15), dots are the simulation 


results elaborated according to Eqs. (16) and (17). 
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Figure 9: (color online) Results of 3D migration defects absorbed by dislocation lines, 
grouped by box sizes with a fixed ratio. (a) sink strength from simulation (empty dots) 
and theoretical expression (dotted and solid lines are the results from Eqs. (18) and (20), 
respectively). (b) the relative error between simulation and theoretical expression (empty 


and solid dots are the results using Eqs. (18) and (20), respectively). 
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Figure 10: (color online) Results of 1D migration defects absorbed by dislocation lines, 
grouped by box sizes with a fixed ratio. (a) sink strength from simulation (empty dots) 
and theoretical expression (solid lines). (b) the relative error between simulation and 
theoretical expression, the solid line represents the error estimated by Eq. (14). (c) data 
obtained by removing the error described by Eq. (14) from (b). Vertical lines are the 


corresponding position of mutation peaks in Fig. 12(b). 
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Figure 11: (color online) Schematics of 1D migration in a periodically extended simulation 
box in OKMC, picture is taken perpendicular to the direction of dislocation lines. Two 
kinds of round circles represent dislocation lines with different radii. The furthest migra- 
tion distance, namely the cutoff distance, is labeled by Deut, and Leut2. Dislocations with 
radius Ra (left) corresponding to a cutoff distance of Leut1, while a slightly reduction of 
the dislocation radius from Rg, to Raz (right) will cause a mutation of the cutoff distance 


from Leuti to Leut: 
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Figure 12: (color online) (a) the mutation of the cutoff distance (in the unit of box size 
lz), Leut, of 1D migration as the function of volume fraction, compared with the average 
migration distance, L, obtained from Eq. (21). (b) the ratio of Leut to L, each peak 
represents a mutation of the cutoff distance. Four high peaks corresponding to the four 
main peaks of relative error in Fig. 10(c), while the four lower peaks corresponding to 


four shoulder peaks of relative error in Fig. 10(c). 
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Figure 13: (color online) The transition from 3D to 1D migrating sink strength of point 
defects absorbed by dislocation lines, grouped by box sizes, where the volume fraction of 
sinks is fixed at f, = 0.005 and fy = 0.03 (denoted by empty and solid dots, respectively). 
(a) simulated sink strength as a function of rotation energy, with temperature T = 723 K. 
(b) the master-curve representation, solid line represents the theoretical relationship of x 
and y in Eq. (15), dots are the simulation results elaborated according to Eqs. (16) and 
(17). 
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